Optical tornadoes in photonic crystals 
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Based on an optical analogy of spintronics, the generation of optical tornadoes is theoretically 
investigated in two-dimensional photonic crystals without space-inversion symmetry. We address 
its close relation to the Berry curvature in crystal momentum space, which represents the non- 
trivial geometric property of a Bloch state. It is shown that the Berry curvature is easily controlled 
by tuning two types of dielectric rods in a honeycomb photonic crystal. Then, Bloch states with 
large Berry curvatures appear as optical tornadoes in real space. The radiation force of such a 
configuration is analyzed, and its possible application is also discussed. 

PACS numbers: 03.65.Vf, 42.15.Eq, 42.70.Qs, 72.25.-b, 



Spintronics is usually discussed in terms of its aspect 
as magnetoelectronics, but it can be also interpreted as a 
technology to control electronic angular momentum via 
band(gap) engineering. The band engineering has been 
a key ingredient in electronics, and nowadays also in op- 
tics [l|, |3| and acoustics [1, |3| , where artificial materials 
with periodic structures play the same role as solid crys- 
tals in electronics. The band engineering in every sys- 
tem realizes various types of functionalities by designing 
bandgaps and band dispersions of Bloch waves. Electron- 
ics utilizes these factors to control the electronic charge 
transport. Then, spintronics attempts to control the elec- 
tronic spin transport. However, in contrast to the charge, 
each spin component is not necessarily conserved in the 
transport process and intricately changes in some cases 
due to the spin-orbit interaction. However, at the same 
time, we can efficiently use this interaction to control 
the spin transport as demonstrated, for example, in the 
spin Hall effect Q , where applied electric field gener- 
ates the spin current perpendicular to the field. Berry's 
geometric phase f?! of a Bloch state is a key ingredient 
for the intrinsic mechanism of this effect [3, 0]. More 
specifically, there appears non-zero Berry curvature in 
each Bloch band relevant to the effect. In most cases, 
the spin-orbit interaction is essential for the emergence 
of this non-trivial geometric property. 

The spin-orbit interaction is a relativistic effect and 
can be crucial not only in a specific crystal structure. 
Photon is a relativistic particle of spin 1 and has the spin- 
orbit interaction in a extreme form which is often rep- 
resented as the transversality condition. Consequently, 
electromagnetic waves in inhomogeneous media also show 
the non-trivial geometric property l3j For instance, 
the transverse shift of an optical beam reflected or re- 
fracted at an interface, which is conventionally known 
as the Imbert-Fedorov effect , in total internal re- 
fiection, is naturally interpreted as the optical spin Hall 
effect The above mechanism of optical Hall 



effect was given a more general flavor on the stage of 
a two-dimensional photonic crystal (PhC) which breaks 
the space-inversion symmet ry ( SIS) [3) [13 and/or the 
time-reversal symmetry [isl . Il9| . When a parameter for 
the symmetry breaking is changed, photonic bands show 
frequent crossing and repulsion with each other. There 
appear large Berry curvatures at nearly degenerate points 
in momentum space [l^]. It was also suggested that a 
wave-packet constructed from Bloch waves around such 
a state seems to involve a kind of rotational motion 17 1 . 



However, any deflnite physical pictures in real space have 
not yet been given to such Bloch waves. 

In this paper, an optical analogy of spintronics in PhCs 
is theoretically investigated. That is, we attempt to uti- 
lize optical angular momentum via band engineering tak- 
ing into account the above background information. For 
this purpose, we address how the Berry curvature in mo- 
mentum space affects physical properties of Bloch waves 
in real space. We focus on the electromagnetic momen- 
tum density and the radiation force of a state with a 
large Berry curvature. It is shown that the Berry curva- 
ture of a photonic Bloch state in momentum space has a 
close relationship to the vortex structure of its momen- 
tum flow in real space. This would open up a new way 
to realize an optical mixer which can selectively scramble 
and transport a specific kind of nanoparticles. It would 
serve for a broad range of research fields as well as the 
optical tweezers technology does. 

Band engineering for the designing of Berry curvature 
— The Berry curvature J7(fc) of a photonic Bloch state 
with a crystal momentum k stands for its geometric as- 
pect [l3] ■ As well as the generality of band engineering 
in several different systems, e.g., electronic, photonic and 
phononic systems, the Berry curvature of a Bloch state 
shares a common property in every system. The most 
prominent one is its scaling behavior around a nearly- 
degenerate point fco ,20]. Especially for the typical case 
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in a two-dimensional system, the scaling behavior is rep- 
resented by 



{v^\k-ko\ 



(1) 



where A is the inter-band separation and w is a nomi- 
nal velocity constant around feg. This property enables 
us to control the Berry curvature via the photonic band 
engineering. 

Here we study this tunability in a honeycomb PhC 
m without SIS. The PhC consists of two (A and B) 
types of circular rods. We can tune the difference be- 
tween the two types of rods through their dielectric con- 
stants £a{b) or radii r^(B). The dielectric constant of 
the background medium Sb is taken to be different from 
both of ea and Eb- When all the rods are identical, the 
PhC recovers the SIS, and its point group is Cqv The 
point group at the K point (Gk) is C31,, which has two 
one-dimensional irreducible representations and one two- 
dimensional irreducible representation. Thus, the modes 
at the K point are either non-degenerate or doubly de- 
generate. When the SIS is broken, the point group of the 
PhC and Gk becomes Gs^ and G3, respectively, and the 
doubly-degenerate modes are lifted. For simplicity, we 
focus on the case ta = tb- The degree of SIS breaking 
and the inter-band separation of the lifted modes at the 
K point are controlled through [sa ^ £b)- 

Berry curvature and optical tornado — Figure [T] shows 
the photonic band structure of a honeycomb PhC with 
the lattice constant a, which is obtained by the photonic 
Korringa-Kohn-Rostoker method The black lines 

stand for the bands of transverse magnetic (TM) mode 
and the red lines for those of transverse electric (TE) 
mode. We can see various types of nearly degenerate 
points. Here we focus on the nearly degenerate partners 
of the TE first and second bands at the K point. The 
Berry curvature of the TE first band is shown in Fig. [51 
As well as the pattern in Fig. ^ the Berry curvature 
of every band has the three-fold symmetry, while every 
band dispersion exhibits six-fold symmetry. For bands 
which are not globally degenerate in the Brillouin zone, 
these symmetric patterns are understood by the following 
arguments. With a symmetry operation S of the point 
group, the energy dispersion a;(fc) satisfies a simple rela- 
tion a;(<S~^fc) = uj{k). On the other hand, the Berry cur- 
vature n(fc) satisfies a relation for an axial vector field 
in momentum space, just like a magnetic field in real 
space. When the time-reversal symmetry is preserved, 
we can also tell uj(—k) = uj{k) and r2(— fe) = — r2(fc). At 
a nearly degenerate point, the Berry curvatures of upper 
and lower bands have the opposite sign with each other 
because of the topological charge (Chern number) con- 
servation [2^. Their signs interchange with each other 
by reversing sa and e b ■ 




FIG. 1: Photonic band structure of the honeycomb PhC with 
ea ~ 1.44, es = l,eb = 4, ta = = 0.2a. Black lines 
for the transverse magnetic (TM) mode and red lines for the 
transverse electric (TE) mode. 
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FIG. 2: Berry curvature of the TE first band in Fig. [U 



Next, let us check the real space configurations of 
states with large Berry curvatures. Figure [3] shows the 
time-averaged Minkowski momentum density D x B oi 
the Bloch state at the K point in the TE first band. 
Only the configuration in a single Wigner-Seitz unit cell 
is shown. The inset is for the state at a point slightly 
shifted from the K point in the same band. In the present 
case, eA{B) < ^b, the momentum flow of the TE first band 
swirls mainly around the center of the unit cell. However, 
in the case eA{B) > ^b, there can appear eddies mainly 
at the corners of the unit cell. In any case, here we call 
such a Bloch state with vortex structure as an optical 
tornado. This is because, just like a real tornado, the 
possible propagation direction of an optical tornado can 
be perpendicular to its rotation vector, as shown in the 
inset of Fig. [3] This point is a noticeable difference be- 
tween the optical tornado and a Laguerre-Gauss beam 
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with orbital angular momentum in free space. 

There is a complementary property between the tor- 
nadoes of upper and lower bands around a nearly degen- 
erate point. The orientation of a tornado around the K 
point is clockwise for the TE first band (Fig. [3]), while 
that is counter-clockwise for the TE second band (not 
shown). Here one may think that the tornadoes are spe- 
cific to the K point in the hexagonal Brillouin zone, which 
is special in a sense that the point group is nontrivial. 
However, we can also find such tornadoes at a nearly- 
accidental-degenerate point of a square lattice PhC with- 
out SIS. The point group of this PhC is trivial, having 
no symmetry operation other than identity one. 
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FIG. 3: Minkowski momentum density of the Bloch state at 
the K point, (|f ,0), of the TE first band, see Fig. [T] Inset 
for the state at the point, (|| — 0.004 x ^,0), of the same 
band. The circles of solid and dashed lines stand for A-type 
and B-type rods, respectively. In each panel, the magnitude 
is normalized by the average energy density/c. Color con- 
tour stands for the absolute value in the main panel and that 
multiplied by 4 in the inset. 



Radiation force induced by optical tornado — Finally, we 
consider the radiation force acting on a small particle in 
the PhC. Here we are assuming the situation in which 
the background medium is an ideal liquid with ef,. When 
the particle size is sufficiently smaller than a typical wave 
length, the particle behaves like an electric dipole Pe in- 
duced by the applied field of a Bloch state. In the non- 
relativistic dynamics of a dipole, the ponderomotive force 
is given in terms of the Lorentz force 2J| as. 



F = {p,.V)E+^xB 

at 



(2) 



The same effect should be also estimated in terms of the 
Maxwell stress tensor T containing both the contribu- 
tions of applied field and dipole radiation as F = J T-ds, 
where the integral is taken over a surface just outside the 
particle. This leads to the estimation in Ref. [25^ which 
is different from Eq. However, in both estimations. 



the time-averaged force takes the common form 26 1 
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where E stands for the complex amplitude of a harmonic 
electric field E{t) = K[£'e-*"*], ap{uj) = aR{uj) + iai{Lu) 
is the complex polarizability of the particle. Readers can 
find the polarizability Q!p(a;) for a sphereical object of a 
radius with a complex permittivity ep(ti-') in Ref. [27| . 
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The last expression of Eq. ([3]) suggests that the net 
force acting on the particle is divided into the gradient 
force (the first term) and the scattering force (the second 
term) . Here it should be noted that the second term also 
includes the force via absorption for particles with inter- 
nal degrees of freedom, e.g., exciton resonance etc. The 
scattering force is further divided into the momentum 
term and the rotational term as follows. 
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where we have used 'V ■ E = 0. It should be noted that 
the rotational term (the last term) vanishes for a simple 
plane wave and TM mode, but not necessarily for TE 
mode. 

For latter convenience, we extend the above argument 
to the case of a line dipole with polarizability per unit 
length ap{uj). The formula for the radiation force (F) can 
be reinterpreted as the radiation force per unit length (F) 
by replacing ap{uj) with ap{u!). For a cylindrical object 
of radius rp with £p{uj), the polarizability per unit length 
is estimated as 



ap{uj) 



- (nr) 
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for TM mode 
for TE mode 



(7) 



,(8) 



where / = 1 for TM mode and Z = 2 for TE mode. The 
denominator of Eq. ([7]) represents the effect of radiation 
damping. 

The above result tells us in which situation the charac- 
teristics of an optical tornado can be efficiently utilized, 
while it should be regarded as an order estimation. For 
the purpose of steady trapping, in which tornadoes un- 
fortunately lose their stage, the gradient force must be 
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dominant. This is the situation required in the optical 
tweezers technology , and is realized under the condi- 
tion aii{u)) ^ cti{uj). Therefore, the particle should not 
show resonant scattering nor absorption at the frequency 
of an applied harmonic field. On the other hand, in or- 
der to utilize the vortex structure of a tornado, we must 
enhance the scattering force. This situation is realized 
under the condition aii{uj) <Si ai^Lu). However, from the 
above formulae, this is not so easy in Rayleigh regime, 
i.e., Tp/X <C 1 where A is the relevant wavelength. We 
need some resonance at the frequency of interest. This is 
because a resonant scattering can help the efficient mo- 
mentum transfer from electromagnetic field to nanopar- 
ticles [2^. In order to investigate this problem more se- 
riously, we pursue a full electromagnetic analysis taking 
into account a resonance. 

To this end, we evaluate the radiation force acting on a 
cylindrical probe via the formula (F) = J (T) ■ ds, where 
the integral is taken over the cylindrical surface of a unit 
length. The Maxwell stress tensor T contains the contri- 
bution of the scattered wave from the cylindrical probe, 
which is calculated by the t-matrix method [sO]. The 
choice of the cylindrical shape for the probe is just for 
a technical simplification in the two-dimensional system. 
In the analysis, a metallic property, i.e., < 0, is intro- 
duced to demonstrate an example of resonant scattering 
due to a mechanism other than exciton excitation. Al- 
though a naive expectation from Eqs. and ^ sug- 
gests the resonance at €p = —4) for TE mode, 
the full electromagnetic analysis for the finite size object 
with Tp = 0.02a predicts the slight shift of the resonant 
point, ep = —4.1. (For a metallic sphere, the resonance 
is expected around Cp = — 2ef,.) Figure |4] shows the re- 
sults, and we can see that the radiation force (per unit 
length) can take over a vortex structure from the optical 
tornado shown in Fig. [3] [3l| . Its physical order of mag- 
nitude is about 2pN/im~^ in the case with a = 1/im and 
the radiation of lmW/im~^. (Here we do not take into 
account the slow- light effect, which enhances the magni- 
tude of the force.) The optical mixer via optical tornado 
is sufficiently feasible when the resonant enhancement of 
the momentum transfer is combined. 

Summary — In summary, we have theoretically investi- 
gated the generation of optical tornadoes and its close 
relation to the Berry curvature in a honeycomb photonic 
crystal without space-inversion symmetry. It was demon- 
strated that the characteristics of an optical tornado is 
easily controlled via band engineering. We also discussed 
the radiation force produced by the tornado acting on 
small objects. An optical mixer via the tornado will be 
realized by combining the resonant enhancement of ef- 
ficient momentum transfer. In other words, this func- 
tionality can be switched selectively for a specific kind of 
nanoparticles. This property would be useful for filtering 
or sorting a mixture of different kinds of particles. 




FIG. 4: Radiation force per unit length induced by the tor- 
nado at the K point of the TE first band in Fig. [T] The probe 
is a thin metallic rod with tp — —4.1 and Vp — 0.02a. The 
magnitude is normalized by the average energy density x a. 
The arrow size is the actual size multiplied by 0.05. 
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